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We consider a power system with N transmission lines whose initial loads (i.e., power flows) 
Li,..., Ln are independent and identically distributed with Pl{x) = P[I/ < a;]. The eapacity Ci 
defines the maximum flow allowed on line i, and is assumed to be given by Ci = (1 + a)Li, with 
a > 0. We study the robustness of this power system against random attacks (or, failures) that 
target a p-fraction of the lines, under a democratie hber bundle-like model. Namely, when a line fails, 
the load it was carrying is redistributed equally among the remaining lines. Our contributions are as 
follows: i) we show analytically that the final breakdown of the system always takes place through 

a first-order transition at the critical attack size p* = 1- /larr ^ i r ^ inn , where E [-1 is 

the expectation operator; ii) we derive conditions on the distribution Pl{x) for which the hrst order 
break down of the system occurs abruptly without any preceding diverging rate of failure; hi) we 
provide a detailed analysis of the robustness of the system under three specihc load distributions: 
Uniform, Pareto, and Weibull, showing that with the minimum load Lmin and mean load E [L] 
hxed, Pareto distribution is the worst (in terms of robustness) among the three, whereas Weibull 
distribution is the best with shape parameter selected relatively large; iv) we provide numerical 
results that confirm our mean-field analysis; and v) we show that p* is maximized when the load 
distribution is a Dirac delta function centered at E[L], i.e., when all lines carry the same load; 
we also show that optimal p* equals This last finding is particularly surprising given that 

heterogeneity is known to lead to high robustness against random failures in many other systems. 

PACS numbers: 64.60.Ht, 62.20.M-, 89.75.-k, 02.50.-r 


I. INTRODUCTION 

As we embark on a future where the demand for elec¬ 
tricity power is greater than ever, and the quality of life 
of the society highly depends on the continuous func¬ 
tioning of power grid, a fundamental question arises as 
to how we can design a power system in a robust and 
reliable manner. A major concern regarding such sys¬ 
tems are the seemingly unexpected large scale failures. 
Although rare, the sheer size of such failures has proven 
to be very costly, at times affecting hundreds of millions 
of people [Il,i ; e.g., the recent blackout in India ii- 
Such events are often attributed to a small initial shock 
getting escalated due to intricate dependencies within a 
power system This phenomenon, also known as 

cascade of failures, has the potential of collapsing an en¬ 
tire power system as well as other infrastructures that 
depend on the power grid [I®; e.g., water, transport, 
communications, etc. Therefore, understanding the dy¬ 
namics of failures in power systems and mitigating the 
potential risks are critical for the successful development 
and evolution of many critical infrastructures. 

In this work, we study the robustness of power sys¬ 
tems under a democratic fiber bundle-like model [llMlSj . 
which is based on the equal redistribution of load upon 
the failure of a power line. It was suggested by Pahwa 
et al. [ll| that equal load redistribution can be a rea¬ 
sonable assumption (in the mean-field sense) due to the 
long-range nature of Kirchoff’s law. This is especially so 
under the DC power flow model that approximates the 
standard AC power flow model when the phase differ¬ 
ences along the branches are small and the bus voltages 
are fixed |ll| . In many cases, power flow calculations 


based on the DC model is known 0111 to give accu¬ 
rate results that match the AC model calculations. 

Our problem setting is as follows: We consider N 
transmission lines whose initial loads (i.e., power flows) 
Li,..., ijv are independently drawn from a distribution 
Pl (x) = P [a < x] . The maximum flow allowed on a line 
i defines its capacity, and is given by Q = {l-\-a)Li with 
a > 0 denoting the tolerance parameter. If a line fails (for 
any reason), its load will be redistributed equally among 
all lines that are alive, meaning that the load carried by a 
line may increase over time. We also assume that any line 
whose load exceeds its capacity will be tripped (i.e., dis¬ 
connected) by means of automatic protective equipments 
so as to avoid costly damages to the system. 

We study the robustness of this system against random 
attacks (or, failures) that target a p-fraction of the lines. 
The failure of the p-fraction of lines may cause further 
failures in the system due to flows of some of the lines 
exceeding their capacity. Subsequently, their load will be 
redistributed which in turn may cause further failures, 
and so on until the cascade of failures stops; note that 
this process is guaranteed to converge, at the very least 
when all lines in the system fail. 

One of our important findings is to show the existence 
of a critical threshold on the attack size p, denoted by p*, 
below which a considerable fraction lines remain func¬ 
tional at the steady state; on the other hand, if p > p*, 
the entire system collapses. We show that the critical at¬ 
tack size is given by p* = 1 - \ l>x])) > 

where E [•] denotes the expectation operator. In addi¬ 
tion, we show that the phase transition at p* is always 
first-order, i.e., the variation of the “fraction of functional 
lines at the steady state” with respect to “attack size p” 
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has a discontinuous first derivative. In a nutshell, what 
this means is that power systems under the democratic 
fiber bundle model tend to exhibit very large changes to 
small variations on the failure size (around p*), render¬ 
ing their robustness unpredictable from previous data. In 
fact, this type of first order phase transition is attributed 
0 to be the origin of large but rare blackouts seen in real 
world, in a way explaining how small initial shocks can 
cascade to collapse large systems that have proven stable 
with respect to similar disturbances in the past. 

Our second main contribution is to demonstrate the 
clear distinction between the case where the first order 
break down of the system occurs abruptly without any 
preceding diverging rate of failure versus the case where 
a second order transition precedes the first-order break¬ 
down. In the former case, ii p < p* the final fraction of 
alive lines will be given by I — p meaning that no single 
additional line fails other than those that are initially at¬ 
tacked, whereas the whole system will suddenly collapse 
if the attack size exceeds p*. These cases are reminis¬ 
cent of the most catastrophic and unexpected large-scale 
collapses observed in the real world. We provide explicit 
conditions on the distribution Pl{x) of the loads and the 
tolerance parameter a that distinguish the two cases. 

Last but not least, we show that p* is maximized when 
the load distribution is a Dirac delta function centered 
at E[L], i.e., when all lines carry the same load. The 
optimal p* is shown to be given by , regardless of the 
mean load E [L]. This finding is particularly surprising 
given that complex networks are known to be extremely 
robust against random failures when their degree distri¬ 
bution is broad [l^ : e.g., when the number of links inci¬ 
dent on a line follows a power-law distribution. 

We believe that our results provide interesting insights 
into the dynamics of cascading failures in power systems. 
In particular, they can help design power systems in a 
more robust manner. The results obtained here may have 
applications in fields other than power systems as well. 
Fiber bundle models have been used in a wide range of 
applications including fatigue [T^ . failure of composite 
materials [l^, landslides [T^, etc. A particularly inter¬ 
esting application is the study of the traffic jams in roads 
[ 23 , where the capacity of a line can be regarded as the 
trafhc flow capacity of a road. 

The paper is structured as follows. In Section El we 
give the details of our system model, discuss how it com¬ 
pares with other models in the literature, and comment 
on its applicability in power systems. Analytical results 
regarding the robustness of the system against random 
attacks are provided in Section Hill for general load distri¬ 
butions. These results are discussed in more details for 
three specific load distributions in Section HVl and various 
load-distribution-specific conclusions are drawn. Section 
El is devoted to numerical results that confirm the main 
findings of the paper for systems of finite size. In Section 
ED we derive the optimal load distribution that leads to 
maximum robustness among all distributions with the 
same mean, and the paper is concluded in Section IVIII 


II. MODEL DEFINITIONS 

We consider a power system with N transmission lines 
whose initial loads (i.e., power flows) Li,...,Ljv are 
independent and identically distributed with Pl{x) := 
P [L < x]. The corresponding probability density func¬ 
tion is given by pl{x) = -^Pl{x). Let Amin denote the 
minimum value L can take; i.e., 

Amin = sup{a; : Pl{x) = 0}. 

We assume that Amin > 0. We also assume that the 
density pl{x) is continuous on its support. 

The capacity of a line defines the maximum power flow 
that it can sustain, and is typically set to be a 

fixed factor of the line’s original load. To that end, we 
let the capacity Ci of line i be given by 

C'i = (l + a)A„ i = (1) 

with a > 0 defining the tolerance parameter. For sim¬ 
plicity, we assume that all lines have the same tolerance 
parameter a, but it would be of interest to extend our 
results to the case where the tolerance parameter ai of 
a line i is randomly selected from a probability distribu¬ 
tion, for each i = 1,..., N. A line fails (i.e., outages) if 
its load exceeds its capacity at any given time. In that 
case, the load it was carrying before the failure is redis¬ 
tributed equally among all remaining lines. 

Our main goal is to study the robustness of this power 
system against random attacks that result with a fail¬ 
ure of a p-fraction of the lines; of course, all the discus¬ 
sion and accompanying results do hold for the robustness 
against random failures as well. The initial set of failures 
leads to redistribution of power flows from the failed lines 
to alive ones (i.e., non-failed lines), so that the load on 
each alive line becomes equal to its initial load plus its 
equal share of the total load of the failed lines. This 
may lead to the failure of some additional lines due to 
the updated flow exceeding their capacity. This process 
may continue recursively, generating a cascade of fail¬ 
ures^ with each failure further increasing the load on the 
alive lines, and may eventually result with the collapse 
of the entire system. Throughout, we let Uooip) denote 
the final fraction of alive lines when a p-fraction of lines 
is randomly attacked. The robustness of a power system 
will be evaluated by the behavior of rioo (p) as the attack 
size p increases, and particularly by the critical attack 
size p* at which noo(p) drops to zero. 

Our formulation is partially inspired by the democratic 
fiber bundle model [l^ , where N parallel fibers with 

random failure thresholds Ci,... ^Cn (i.e., capacities) 
drawn independently from Pc ( 5 ) share equally an ap¬ 
plied total force of F ; see also [^, [25l - lT^. This model 
has been recently adopted by Pahwa et al. in the con¬ 
text of power systems with F corresponding to the total 
load that N power lines share equally. A major difference 
of our setting with the original democratic fiber-bundle 
model is that in the latter the total load of the system 
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is always fixed at F. This ensures that the load that 
each alive line carries at any given time is independent of 
the specific set of lines that have failed until that time. 
For example if M lines out of the original N are alive, 
one can easily compute the load per alive line as F/M re¬ 
gardless of which N — M lines have actually failed. In our 
model, however, the initial loads of N lines are random 
and they differ from each other, and so do their capaci¬ 
ties. This leads to strong dependencies between the load 
of an alive line and the particular N — M set of lines 
that have failed, and makes it impossible to compute the 
former merely from the number of failed lines. For in¬ 
stance, at any given time, lines that are alive are likely 
to have a larger capacity, and thus a larger initial load 
in view of (HD, than those that have failed. In addition, 
the total load shed on to the alive lines is not given by 
(iV —M)E [L], since the lines that have failed are likely to 
have a smaller capacity, and thus a smaller initial load, 
than average. As a result of these intricate dependencies, 
analysis of cascading failures in our setting becomes sub¬ 
stantially more challenging than that in the fiber-bundle 
model; see Section nil for details. 

We believe that our problem formulation can lead to 
significant insights for the robustness of power systems 
(and possibly of other real-world systems) that can not 
be seen in the original fiber-bundle model. First of all, 
our formulation allows analyzing the robustness of the 
system against external attacks or random line failures, 
which are known to be the source of system-wide black¬ 
outs in many interdependent systems ; the 

standard fiber-bundle model is instead concerned with 
failures triggered by increasing the total force (i.e., load) 
applied to the system. Secondly, unlike the democratic 
fiber bundle model where all lines start with the same 
initial load [1^, power lines in real systems are likely to 
have different loads at the initial set-up although they 
may participate equally in taking over the load of those 
lines that have failed; intuitively speaking, this is also the 
case for traffic flow on roads. 

Our model has some similarities also with the CAS¬ 
CADE model introduced by Dobson et al. [1^. There, 
they assume that initial loads Li,..., are uniformly 
distributed over an interval (imin,Amax), and all lines 
have the same capacity C = This is a signihcant 

difference from our model where capacities vary accord¬ 
ing to (HD. Another major difference is that in the CAS¬ 
CADE model, a fixed amount A is redistributed to all 
alive lines irrespective of the load being carried before 
failure. Therefore, strong dependencies between particu¬ 
lar lines failed and the load carried by alive lines do not 
exist in the CASCADE model. 

A word on notation in use: The random variables (rvs) 
under consideration are all defined on the same probabil¬ 
ity space (D,J^, P). Probabilistic statements are made 
with respect to this probability measure P, and we de¬ 
note the corresponding expectation operator by E. The 
indicator function of an event A is denoted by 1 [A]. 


III. ANALYTIC RESULTS 
A. Recursive Relations 


We now provide the mean-field analysis of the cas¬ 
cading failures of lines for the model described in Sec¬ 
tion El We start by deriving recursive relations concern¬ 
ing the fraction ft of lines that are failed at time stage 
t = 0,1,.... The number of links that are still alive at 
time t is then given hy Nt = N{1 — ft) for alH = 0,1,.... 
The cascading failures start with a random attack that 
targets a fraction p of power lines, whence we have /o = p. 
Upon the failure of these fop lines, their load will be redis¬ 
tributed to the remaining {1 — fo)N lines. The resulting 
extra load per alive line, Qo is given by 


Qo 


E[L]pN 

(l-p)iV 


E|L1^ 


( 2 ) 


At this initial stage, since the pN lines that have been 
attacked are selected uniformly at random, the mean to¬ 
tal load that will be transferred to the remaining lines is 
just given by E [L] pN. 

Now, in the next stage a line i that survived the initial 
attack will fail if and only if its new load reaches its 
capacity [3l|; i.e., if 


Li -l- Qo ^ (1 + Oi)Li, 


or, equivalently if Li < Qo/a. Therefore, at stage t = 1, 
an additional fraction P [A < Qo /ct] of lines will fail from 
the lines that were alive at the end of stage 0. This gives 


fi — /o + (l —/o)P [L < Qo/a] 


l-(l-/o)P 



In order to compute Qi, i.e., the total extra load per 
alive line at stage 1, we should sum the total load of all 
the failed lines until this stage and divide it by the new 
system size 1 — /i. So, Qi(l — /i) is given by the sum of 
Qo(l — /o), and the total load of the lines failed at stage 
1 normalized by the number of lines N' i.e., of lines that 
survived the initial attack but have load L < Qola. Let 
A be the initial set of lines attacked. We get 

i^A-.LiKQoja 

J2Lil[L,<Qo/a] 
i^A 

= Qo(l - /o) + ^ ■ X! ® ^ 

i^A 

= Qo(l - fo) + (1 - /o)E [Al [A < Qo/a]], 
where the last step uses |A|/A = p = fo. Thus, we get 
^ _pE [A] + (1 - p)E [A ■ 1 [A < Qo/a]] 


Qi(l - /i) — Qo(l - /o) + ^ ■ E 
= Qo(l-/o) + ^-E 











4 


upon noting ([2]). We find it useful to note that 

E [L • 1 [L < Qo/a]] = E [L | L < Q^/a] P [L < Q^/a]. 

The general form of ft and Qt will become apparent as 
we compute them at stage t = 2. This time we argue as 
follows. For a line to still stay alive at this stage, two con¬ 
ditions need to be satisfied: i) it should not have failed 
until this stage, which happens with probability 1 — /i; 
and ii) its load should satisfy L > Qija so that its ca¬ 
pacity is still larger than its current load. One additional 
note is that a line that satisfies condition (i) necessarily 
have a load L > Qo/a. Collecting, we obtain 

/2 = l-(l-/i)P[T>Qi/a I T>go/a]. 

The total load that will be redistributed to the remaining 
lines can then be computed as before: 

Q 2 (l-/ 2 )=go(l-/o) + ^-E ^ U 

i^A\Li<Q\ j a. 

= 0o(l - /o) + (1 - /o)E [LI [L < Qi/a]] . 

One can complicate the matters a little bit and get the 
same expression by writing 


e2(i-/2) = 0i(i-/i) + 4E 


i^A-.Q^ ! OL<Li<Q\ j ol 


as well. 

The form of the recursive equations is now clear. Let 
/o = p, No = N{1 - p), and Qo =E[L]jf-. For conve- 


have 


L > ^ 


ft+i= l-(l-/t)P 

pElL] + {l-p)E[L-l[L<^] 

Qt+i= - 137 ^ 7 ^-^ 

Nt+i = (1 - ft+i)N 


L > 


Qt-1 


(3) 


From © we see that cascades stop and a steady is 
reached, i.e., iVt +2 = Nt+i, if 


\l > 

L> ^ 

a 

a 


= 1 . 


(4) 


In order to understand the conditions that would lead to 
©2 we need to simplify the recursion on ft . This step is 
taken in the next section. 


Applying these recursively, we obtain 

t 

1 - ft+i = {1 - fo)l[F[L > Qe/a \ L > Q,_i/a], 
£=0 

where Q-i = 0 as before. Since Qt is monotone increas¬ 
ing in t, i.e., Qt+i > Qt for all t, we further obtain 

1 - ft+i 

= (l-/o)- 


L> ^ 


L > 

Qo 


Qt-1 


p 

L > 

06 

P 

L > ^ 

06 

p 

L > 

06 

P 

L > So 

06 


L > 


= {l-fo)E[L>Qt/a] 


(5) 

Reporting this into © and recalling that /o = p, we 
get the following simplified recursions: 


L > % I L > 


ft+i = 1 - (1 - /t 

^ _ pEIL] + {1-p)E\l\[L<^]i 

^*+1 “ (l-p)P[L>Qt/a] 

Nt+i= {l-p)¥[L>Qt/a]N 


Qt-1 


( 6 ) 


we 

F 

pE[L]^{l-p)E 

T \ 

LI 

L < Si 
- 06 





a(I - p)P [L > Qt/a] 

a 


Failures will stop and a steady-state will be reached when 
ft +2 = ft+i- From the first relation in dH]), we see that 
this holds if 

¥[L>Qt+i/a\ L>Qt/a] = l, 
or, equivalently if 


= 1, (7) 

as we use the middle equation in (jb]). 

Define x := Qtjct, and realize that 

pE[L] + {1 - p)E[L ■ 1[L < x]] 

= pE [L] + (1 - p)E [L • (1 - 1 [L > *])] 

= E[L]-(l-p)E [L-l[L>a:]]. 

With these in place, the condition for cascades to stop 
© gives 


L > 


E[L]-(l-p)E [L-1[L>3;]] 
a(l — p)P [L > x] 


L > X 


= 1 . 


( 8 ) 


B. Conditions for steady-state via a simplification 

Applying the first relation in ([3|) repeatedly, we see 
that 

1 - ft+i = (1 - ft)E [L > Qt/a I L > Qt-i/a] 

^ — ft = (1 ~ ft-iW > Qt-ilct I L > Qt-^l o] 

l-/i =(l-/o)P[L>Qo/a] 


It is now clear how to obtain the final fraction of power 
lines that are still alive at the end of the cascading fail¬ 
ures: One must find the smallest solution X* of ©. Then, 
the final fraction noo{p) of alive lines is given (see ©) by 

riooip) = 1 - /oo = (1 - p)E [L > X*] . (9) 

Under the enforced assumptions on the distribution of 
L, we see that ([5]) holds in either one of the following 
cases: 
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If * > K[L]-(l-p)EfL.l[L>x]] . 

— a{l-p)P\L>x\ > 

ii) If r < nL]-{i-pW[L-l[L>x]] , 

ll) il a; < a{l-p)¥[L>x] 

E [L] — (1 — p)E [L ■ 1[L > x]] 


L > 


a(l — p)P [L > x] 


= 1 . ( 10 ) 


We see that in the latter case, it automatically holds 
P [i > a;] = 1, meaning that the final system size equals 
1 — p. In other words, no single line fails other than the 
pN lines that went down as a result of the initial attack. 
Using P [U > a;] = 1 in (fTUl) . we see that this happens 


whenever 


L > 


pnL] 


a(l-p) 


= 1, which can be regarded as 


the condition for no cascade of failures. This condition 
can help in capacity provisioning, i.e., in determining the 
factor a needed for robustness against p-size attacks, and 
can be rewritten as 



(a) 


> 


pE[L] 
a(l -p)' 


( 11 ) 


The first condition, on the other hand, amounts to 


P [L > x] (ax + E [T 


IN ]E[L] 
L>x])> 


P 


( 12 ) 


We can now see that the final system size noo{p) is al¬ 
ways given by (1 — p)P [L > x*] where x* is the small¬ 
est solution of (ini). This is clearly true for the case (i) 
given above. To see why this approach also works for 
the case (ii), observe that when dm) holds (na) is sat¬ 
isfied for any x in ; .^min]- Hence, the smallest 

solution X* of ca will always give x* < Tmin, leading to 
(1—p)P [L > X*] = 1—p. As discussed before, no cascade 
takes place under m (i.e., in the case {ii) above), so the 
final system size is indeed 1 — p. 

For a graphical solution of noo(p), one shall plot 
P [L > x] (ax -I- E [L I L > x]) as a function of x (e.g., see 
Figure 1(a)), and draw a horizontal line at the height 
E [L] /(I — p) on the same plot. The leftmost intersec¬ 
tion of these two lines gives the operating point x*, from 
which we can compute noo(p) = (1 — p)P [L > x*]. When 
there is no intersection, we set x* = oo and understand 
that noo{p) = 0. 


C. Rupture Condition 

We now know how to compute the final system size 
noo{p) for a given attack size p. In many cases, we will 
be interested in the variation of noo(p) as a function of p. 
This will help us understand the response of the system 
to attacks of varying magnitude. Of particular interest 
will be to derive the critical attack size p* such that for 
any attack with size p > p*, the system undergoes a 
complete breakdown leading to noo(p) = 0. 



(b) 


FIG. 1. (Color online) We demonstrate the distinction be¬ 
tween an abrupt Grst-order rupture, and a Grst-order rupture 
that is preceded by a diverging failure rate. Pl{x) is assumed 
to be of uniform density over the range [Lmin, Tmax] = [10, 50]. 
In both plots, Red curves stand for the case where a — 0.2, 
whereas Blue curves represent a = 1.2. Fig ure U(a)| shows 
V[L > x] {ax -I- E [L | L > x\), whereas Figure 1(b) plots the 
corresponding variation of Uoofp) with attack size p. We ob¬ 
serve that for a = 0.2 (Red), P [L > x] {ax -f E [L | L > x]) 
takes its maximum at the point x = Lmm = 10. As a result, 
we see an abrupt first-order transition of n^o (p) as it suddenly 
drops to zero at the point p = p* = 0.0625, while decaying 
linearly as 1 —p up until that point. The case where a = 1.2 is 
clearly different as P [L > x] {ax -|- E [L | L > x]) is now max¬ 
imized at X = 17.6 > I/min- As expected from our discussion, 
this ensures that the total failure of the system occurs after 
a diverging failure rate is observed. This divergence is clearly 
seen in Figure p (b)| where the dashed line corresponds to the 
1 — p curve. 
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From (IT^ and the discussion that follows, we see that 
the maximum attack size p* is related to the global max¬ 
imum of the function P [L > a:] {ax +K[L \ L > x]). In 
fact, it is easy to see that 


p* = l- 


E[L] 


max{P [L > x] {ax -I- E [L | L > a;])} 


(13) 


The critical point x* that maximizes the function 
P [L > a;] (aa; -I- E [L | L > x]) can shed light on the type 
of the transition that the system undergoes as the at¬ 
tack size increases. First of all, the system will al¬ 
ways undergo a first-order (i.e., discontinuous) transi¬ 
tion at the point p*. This can be seen as follows: We 
have noc{p* ) = 0 by virtue of the fact that no x will 
satisfy (HU, and cascading failures will continue until 
the whole system breaks down. On the other hand, 
noo{p* ) = (1 ~ p)P > 2 :*] > 0 where x* is the point 
that maximizes P [L > a;] {ax -I- E [L | L > x]). We can 
see why it must hold P [L > a;*] > 0 via contradic¬ 
tion: P [T > a;*] = 0 implies that the maximum value of 
P [L > a;] (aa: -I- E [L | T > a:]) is zero, which clearly does 
not hold since at a; = 0 this function equals E [L] > 0 by 
non-negativity of L. 

An interesting question is whether this first order rup¬ 
ture at the point p* will have any early indicators at 
smaller attack sizes; e.g., a diverging failure rate leading 
to a non-linear decrease in noo{p)- With Lmin > 0, we 
know from m that for p sufficiently small, there will be 
no cascades and noo(p) will decrease linearly as 1—p. This 
corresponds to the situations where (fT^ is satisfied at a 
point X < Lmin, i-S., when P [L > a;] {ax -I- E [T | L > x]) 
is linearly increasing with x. An abrupt first-order tran¬ 
sition is said to take place if the linear decay of rioo (p) is 
followed by a sudden discontinuous jump to zero at the 
point p*. Those cases are reminiscent of the real-world 
phenomena of unexpected large-scale system collapses; 
i.e., cases where seemingly identical attacks/failures lead¬ 
ing to entirely different consequences. 

It is easy to see that an abrupt transition occurs if 
P [L > a;] {ax -I- E [L | L > x]) takes its maximum at the 
point x = Amin; §60 Figure [T] In that case, (IT^ either 
has a solution at some x < Amin so that noo{p) = I — P, 
or has no solution leading to naa{p) = 0. Under the as¬ 
sumptions enforced here, P [A > a;] {ax -I- E [A | L > x\) 
is continuous at every a; > 0. Given that this function 
is linear increasing on the range 0 < a; < Amim a maxi¬ 
mum takes place at a; = Amin if at that point the derivate 
changes its sign. We have 

-^(P[A>a:](aa:-bE[A I A>x])) (14) 

dx 

= (axP [A > x] -I- E [A • 1 [A > x]]) 
dx 

= aP [A > x] -I- ax(— pl ( x )) ^ (^J tPL{t)dt^ 

= aP [A > x] -I- ax{—pL{x)) — xpl{x) 

= aP [A > x] — xpL{x){a -|- 1) (15) 


where in the second to last step we used the Leibniz 
integral rule. As expected, for x < Amin, we have 
P [A > x] = 1 and pl{x) = 0, so that the derivative is 
constant at a. For an abrupt rupture to take place, the 
derivative should be negative at the point x = Amin; i-e., 
we need 


a Amin ■ Pl/(Amin)(Q; -f 1) < 0, 
or, equivalently 


, , I \ r < Pi (Amin). (16) 

{a IjAmin 

It is important to note that dm) ensures only 
the existence of a local maximum of the function 
P [A > x] (ax -I- E [A I A > x]) at the point x = Amin- 
This in turn implies that there will be a first order jump 
in Uoo(p) at the point where E [A] /(I—p) = aAmin+E [A]; 
i.e., at the point p that satisfies (jTT]) with equality. How¬ 
ever, for this condition to lead to an “abrupt” first-order 
breakdown, we need x = Amin to be the global maxi¬ 
mum. This can be checked by finding all x that make 
the derivate at (fTSl) zero, and then comparing the corre¬ 
sponding maximum points. If x = Amin is only a local 
maximum, then the system will have a sudden drop in 
size at the corresponding attack size, but will not un¬ 
dergo a complete failure; the complete failure and the 
drop of noo(p) to zero will take place at a larger attack 
size where, again there will be a first-order transition; 
e.g., see Figure^] 

We close by giving the general condition for first-order 
jumps to take place. We need a change of sign of the 
derivative at im, leading to 


±aP [A > x] 


xpL(x)(a -I- I) 


X — X 


.± 


< 0 . 


Equivalently, a first-order jump will be seen for every x* 
satisfying 


Pl{x* ) < 


aP [A > X*] 
(a -b l)x* 


< Pl{x*^). 


(17) 


IV. RESULTS WITH SPECIFIC 
DISTRIBUTIONS 

We analyze a few specific distributions in more details. 
Namely, we will consider Uniform, Pareto, and Weibull 
distributions. 


A. Uniform distribution 

Assume that loads Ai,...,A^r are uniformly dis¬ 
tributed over [Amin, Amax]- In other words, we have 

Pl{x) = — r * 1 [Amin ^ X Amax] , 

-^max -^min 
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so that 

¥[L>x]= 1 < a; < L„,ax] 

-^max -^min 

+ 1 [a; < Linin] • (18) 

We see that over the range x in [0,iniax), the derivative 
of P [L > a;] (aa; + E [L | L > x\) (see (fT^ i is either never 
zero or becomes zero only once at 

For the latter to be possible, we need 2 a+i -^max > 

If the opposite condition holds, i.e., if 2 (^i -^max < Fmin, 
then P [F > a:] (aa; + E [L | L > x\) is maximized at x = 
-bmin, and an abrupt first order break down will occur 
(as p increases) without any preceding diverging failure 
rate. As expected, the condition 2aJr\ ^™ax < L-^ain is 
equivalent to the general rupture condition (1161) and can 
be written most compactly as 

Fniin 

a < - - - 

max (Aniax 2Z/niin, 0) 

It follows that if Amax < 2Lmin, then an abrupt rupture 
takes place irrespective of the tolerance factor a. 


B. Pareto distribution 


Distribution of many real world variables are shown 
to exhibit a power-law behavior, with very large variabil¬ 
ity [s^ - lssjl . To consider power systems where the initial 
loads of the lines can exhibit high variance, we consider 
the case where Li,..., are drawn from a Pareto dis¬ 
tribution: Namely, with 6, Amin > 0, we set 

PL^^') — 1 [ 3 : ^ • 


To ensure that E \L] is finite, we also enforce that 6 > 1; 
in that case we have E [L] = . Then, the condition 

for an abrupt first order rupture (1161) gives 


(a -|- l)Tn 


r b 


- -b-1 
'^min ’ 


or, equivalently < b. With 6 > 1, this always holds 
meaning that when the loads are Pareto distributed, 
there will always be an abrupt first order rupture at the 
attack size »* = 1 — — = 1-In fact, 

r E[L]-|-ctLmin i+aS-jJ- ’ 

we can see that this attack will lead to a complete break¬ 
down of the system since for x > Amin; we have 


— (P [L > x] {ax -I- E [T I L > x])) (19) 

= aP [L > x] — xpL{x){a + 1) 

= - (a + l)xL^„6x"''-^ 

= {a - b{a + 1)) 

< 0 


for any a > 0 and 6 > 1. Therefore, it is always the 
case that (P [L > x] (ax -I- E [L | L > x])) has a unique 
maximum at x = Amin, and the abrupt first order rupture 
completely breaks down the system. 

These results show that for a given Lmin and E [L] with 
E [L] > Lmin, Pareto distribution is the worst possible 
scenario in terms of the overall robustness of the power 
system. Put differently, with Lmin and E [L] fixed, the 
robustness curve nao{p) for the Pareto distribution con¬ 
stitutes a lower bound for that of any other distribu¬ 
tion. From a design perspective, we see that changing 
the tolerance parameter a will not help in mitigating the 
abruptness of the breakdown of the system in the case of 
Pareto distributed loads. On the other hand, the point 
at which the abrupt failure takes place, i.e., the critical 
attack size p* can be increased by increasing a. 


C. Weibull distribution 


The last distribution we will consider is Weibull distri¬ 
bution, which has the form 


Pl{x) = ^ 


-^min 


^ Tniin] 


with A, fc > 0. The case k = 1 corresponds to 
the exponential distribution, and k = 2 corresponds 
to Rayleigh distribution. The mean load is given 
by E [L] = Lniin + Ar(I + 1/k), where r(-) is the 
gamma-function. As usual, we check the derivative of 
(P [L > x] (ax + E [L I L > x])) for x > Amin- On that 


^min 


so that 


dx 


(P [L > x] (ax + E [A I A > x])) 


= e 


a — (a + I)x 


k f X — L-a 


A 


A 


fc-i' 


( 20 ) 


which becomes zero if 


x(x Amin) " — 




(a -I- l)k 


( 21 ) 


This already prompts us to consider the cases k < 1 
and fc > 1 separately. In fact, with /c > I, we see that 
AL(Amin) = 0 and (fTOj) does not hold regardless of a. In 
addition, there is one and only one x > Amin that can 
satisfy (ED. Consequently, for k > 1 the system will al¬ 
ways undergo a second-order transition with a diverging 
rate of failure before breaking down completely through 
a first-order transition. 

The case k < 1 gives an entirely different picture since 
PL(Amin) = oo and (ITSl) always holds regardless of a. 
So, the system will always go through an abrupt first 
order transition at the attack size p* = 1 — —. 

Whether this rupture will entirely breakdown the system 
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FIG. 2. The two stage breakdown of the system is demon¬ 
strated, where Li,..., Ln are drawn from Weibull distribu¬ 
tion with k = 0.8, A = 150, Lmin = 10, a = 0.2. We plot 
the relative final size noo(p) as a function of the attack size 
p. The Inset zooms in to the region where the system goes 
through a series of first-order, second-order, and then again a 
first-order transition. 


depends on the existence of the solutions of (1^ . It is 
easy to see that x(x — takes its minimum value 

at a; = Lmin/k and equals to ^ ■ Thus, if it 

holds that 


k ) a + 1 


( 22 ) 


then m has no solution and the derivative given 
at dini) is negative for all x > meaning that 

P [L > a;] (ax + E [L | L > a;]) is maximized at a; = Amin- 
Then, the abrupt hrst order rupture at p* = 1 — 
E[L]+aL indeed breakdown the system completely. 

The same conclusion follows if (l2^ holds with equality 
by virtue of the fact that (EUl) is again non-positive for 
all X ^ Amin- 

On the other hand, if 


k ) ^ a -f 1 


(23) 


then m will have two solutions both with x > Amin- 
This implies that P [A > a:] (era; -I- E [A | A > a;]) has an¬ 
other maximum at a point x > Amin- If this maximum 
is indeed the global maximum (i-e., it is larger than the 
maximum attained at a; = Amin), then the system will 
go under two first-order phase transitions before break¬ 
ing down. First, an abrupt rupture will take place at 
p* = 1 — ■ But, this won’t break down the 

system completely and noo(p*^) will be positive. As p 
increases further, we will observe a second-order transi- 


FIG. 3. (Color Online). We plot noo(p) vs. p under six 
different cases. Analytical results are represented by lines, 
whereas empirical results (obtained through averaging over 
500 independent runs) are represented by symbols. We set 
N — 100,000, Aniin = 10, and E [A] = 30. For the case 
when Ai,..., Ljv follow a Weibull distribution, we take the 
shape parameter to be k — 2, leading to a scale parameter 
A = 22.5676. We see that numerical results match the ana¬ 
lytical results very well. 


tion with a diverging rate of failure until another first- 
order rupture breaks down the system completely. We 
demonstrate this phenomenon in Figure [21 where we set 
k = 0.8, A = 150, Amin = 10, a = 0.2. We empha¬ 
size that this behavior (i.e., occurrence of two first-order 
transitions) is not immediately warranted under (1231) . It 
is also needed that P [A > a;] (as -|- E [A | A > s]) has a 
global maximum at a point x > Amin. 


V. NUMERICAL RESULTS 

We now check the validity of our mean-field analysis 
for finite number N of power lines via simulations. We 
will do so with an eye towards comparing the robustness 
of power systems under different distributions of loads. 

In the first batch of simulations, we fix the minimum 
load at Amin = 10 and mean load at E [A] = 30. These 
constraints fully determine the load distribution pl(x) 
in the cases where pl is Uniform (with Amin = 10 and 
Amax = 50) or Pareto (with Amin = 10 and b = 1.5). For 
the case where pp is Weibull, we need to pick k and A 
such that Ar(l -b 1/fc) = 20, where r(a;) = 

We consider k = 2,X = 22.5676 as an example point. 

Our simulation set up is as follows. We fix the 
number of lines N, and generate N random variables 
from the given distribution pl{x) corresponding to loads 
Ai,..., Ljy. Then, for a givenp, we perform an attack on 
[piV] lines that are selected uniformly at random and as- 
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sume those lines have failed. Next, using the democratic 
redistribution of loads, we iteratively fail any line whose 
load exceeds its capacity (which is set to (1 + 0 ) times its 
initial load). We consider two possible tolerance param¬ 
eters: i) a = 0.2 and ii) a = 0.7. The process stops when 
the system is stable; i.e., all lines have a load below their 
capacity. Of course, this steady state can be reached at 
a point where all lines in the system have failed. We 
record the corresponding fraction of active lines at the 
steady state. This process is repeated independently 500 
times for each p, and the average fraction of active lines 
at the steady state over 500 independent runs gives the 
empirical value of noo{p)- We then compare this with the 
quantity obtained from our analysis in Section IlIIl 

Results are depicted in Figure O First of all, we see an 
almost perfect agreement between our mean-held analy¬ 
sis and numerical results. It is worth noting that the Ht 
between analysis and simulations required particularly 
large values of N in the case of Pareto distribution. For 
the other two distributions, even N = 5000 leads to al¬ 
most perfect agreement. Focusing on two curves corre¬ 
sponding to uniform distribution, we see from Figure [3] 
that the tolerance parameter a not only changes the max¬ 
imum attack size p* that the system can sustain (in the 
sense of not breaking down entirely), but it can also af¬ 
fect the type of the phase transition. In particular, with 
a = 0.2, an abrupt failure takes place at p* = 0.0625, 
whereas with a = 0.7 the system goes through a second 
order transition starting with the attack size p = 0.189, 
and then breaks down entirely through a first-order jump 
at p* = 0.203. 

As expected from our previous discussion, the distribu¬ 
tion that leads to the worst robustness is Pareto among 
all distributions considered here. However, we see that 
under certain conditions uniform and Weibull distribu¬ 
tions can match the poor robustness characteristics of 
the Pareto distribution; one example is the case shown 
in Figure [ 3 ] with uniform load distribution and a = 0.2. 
Finally, we observe that Weibull distribution can lead to 
a significantly better robustness than Pareto and Uni¬ 
form distribution, under the same mean and minimum 
load. 

The last observation worths investigating further. In 
particular, even when Lmin and E [L] are fixed, the 
Weibull distribution has another degree of freedom; i.e., 
parameters k and A are arbitrary subject to the condition 
that Ar(l -I- 1/k) = E [L]. In order to understand the ef¬ 
fect of the shape parameter k in the robustness of power 
systems under Weibull distributed loads, we ran another 
set of simulations with N = 100, 000, a = 0.7, Amin = 10, 
E [L] = 30, and A = various values of k. 

Results are depicted in Figure |4] where again analytical 
results are represented by lines and empirical results (ob¬ 
tained through averaging over 500 independent runs) are 
represented by symbols. We again observe an excellent 
match between analytical and numerical results. We re¬ 
mark that with the given parameter setting, the cases 
where A: < 1 all result in the same robustness behavior 



FIG. 4. (Color Online). We plot noo(p) vs. p when Li,..., Ljv 
follow a Weibull distribution with Lmm = 10, E [L] = 30. We 
set N = 100, 000 and a = 0.7. Analytical results are repre¬ 
sented by lines, whereas empirical results (obtained through 
averaging over 500 independent runs) are represented by sym¬ 
bols. Again, we see that numerical results match the analyt¬ 
ical results pretty well. 


with an abrupt first-order rupture at p = 0.189. 

More importantly, we see that the robustness of the 
system improves as the parameter k increases. It is 
known that as k gets larger the Weibull distribution gets 
closer and closer to a Dirac delta distribution centered at 
its mean. In other words, as k goes to infinity the Weibull 
distribution converges to a degenerate distribution and 
loads Ai,...,Ajv will all be equal to the mean E[L]. 
This naturally prompts us to ask whether a degenerate 
distribution of loads is the universally optimum strat¬ 
egy among all possible distributions with the same mean 
E [L] , with optimality criterion being the maximization 
of robustness against random attacks or failures. Here, a 
natural condition for maximization of robustness would 
be to maximize the critical attack size p*. We answer 
this question, in the affirmative, in the next section. 


VI. OPTIMAL LOAD DISTRIBUTION 

To drive the above point further and to better under¬ 
stand the impact of the shape parameter k on the sys¬ 
tem robustness, we now plot the maximum attack size 
p* as a function of k under the same setting; see Figure 
[S] Namely, we let Li,... ,Ln follow a Weibull distribu¬ 
tion with Amin = 10, and A = 20/F(l -|- 1/k) so that 
E [A] = 30. We see from Figure [S] that, in all choices of 
a considered here, the maximum attack size p* is mono¬ 
tone increasing with fc; note that p* is seen to be constant 
over the range 0 < k < 1. It is also evident from Figure 
[S]that p* tends to converge to a fixed value as fc —>■ 00 . 















. attack size, 
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FIG. 5. (Color Online). We plot the maximum attack size p*, 
when Li,..., Ln follow a Weibull distribution with Lmin = 
10, E [L] — 30, as a function of the shape parameter k of 
the Weibull distribution. We set N = 100, 000 and consider 
four tolerance parameters a = 0.1, 0.3, 0.5, 0.7. The curves 
correspond to analytical results computed directly from 03). 


On the other hand, with /c —>■ oo, we know that Weibull 
distribution converges to a Dirac delta distribution cen¬ 
tered at E [L]. It is therefore of interest to check whether 
p* is always maximized by choosing all loads Li,..., Ljv 
equally, i.e., by choosing pl{x) to be a degenerate distri¬ 
bution with mean E [L] and zero variance. 

Let pl{x) be an arbitrary distribution with mean E [L], 
and assume that pl(x) = 0 for cc < 0; i.e., that L is non¬ 
negative. Recall that maximum attack size p* is given 
by (fT3) and observe that 

F[L > x] (ax -I- E [L | L > x\) 

= axF [L > x] +F[L ■ 1[L > a;]] 
<aE[L]+E[L-l[L> a;]] (24) 

<{a + 1)E [L], (25) 

for any a; > 0. In (1241) we used the Markov Inequality [H, 
p. 151], i.e., the fact that P [L > a;] < E [L] jx for any 
non-negative random variable L and a; > 0. Reporting 
(I25() into m, we get 


p* < 1 - 


E[L] 


(a + 1)E [L] a -I- 1 


(26) 


This shows that the maximum attack size can never 
exceed under any choice of load distribution. On 
the other hand, consider the case where pl{x) = (5(E [L]) 
with (5(-) denoting a Dirac delta function. This implies 
that Li = ■ ■ ■ = L]y = E[L]. Let denote the cor¬ 

responding maximum attack size. With x = E [L] , we 
have P [L > x] = 1 and E [L • 1 [L > xj] = E [L\. Thus, 


lim axP [L > x] -|- E [L ■ 1 [L > x]] = (a -|- 1)E [L] 

x'[W.[L] 


SO that 


max{P [L > x] (ax -I- E [L | L > x])} > (a -I- 1)E [L]. 

X 


Invoking (fT3) . this leads 


I^dirac — 


a 

a -I- 1 


But, (l26ll holds for any distribution and hence is also 
valid for p^irac- Combining these, we obtain that 


Pdh 


a 

a -|- 1 


(27) 


This establishes that a degenerate distribution is indeed 
optimal for any given mean value of the load, and the 
achieved maximum attack size is given by a/(a + 1). 
What is even more remarkable is that, this maximum 
attack size is independent of the mean load E [L]. 

It is now clear to what point the curves in Figure [Sj 
tend to converge as fc —>■ oo; they can indeed be seen 
to get closer and closer to the corresponding value of 
a/(a + 1). We close by demonstrating the variation of 
the final system size as a function of the attack size, in the 
case where loads follow a Dirac distribution. We easily 
see that P [L > x] (ax -I- E [L | L > x]) increases linearly 
for X < E [L] and equals to zero for x > E [L]. Therefore, 
the breakdown of the system will always be through an 
abrupt first order rupture. 

This is demonstrated in Figure [51 where it is seen once 
again that numerical results match the analysis perfectly. 
Comparing these plots with Figures [3] and 13 we see the 
dramatic impact that the load distribution has on the ro¬ 
bustness of a power system. For instance, with a = 0.2 
and mean load fixed at 30 we see that maximum attack 
size that the system can sustain is 6.3% for Pareto and 
Uniform distributions whereas it is 17% when all loads 
are equal. Similarly, with a = 0.7 we see that maximum 
attack size is 18% for Pareto distribution and 19% for 
Uniform distribution, while for the Dirac delta distribu¬ 
tion, it increases to 41%. These findings suggest that 
under the democratic fiber bundle-like model considered 
here, power systems with homogenous loads are signifi¬ 
cantly more robust against random attacks and failures, 
as compared to systems with heterogeneous load distri¬ 
bution. 


VII. CONCLUSION 

We studied the robustness of power systems consisting 
of N lines under a democratic-fiber-bundle like model 
and against random attacks. We show that the system 
goes under a total breakdown through a first-order tran¬ 
sition as the attack size reaches a critical value. We de¬ 
rive the conditions under which the first-order rupture 
occurs abruptly without any preceding divergence of the 
failure rate; those situations correspond to cases where 
no cascade of failures occurs until a critical attack size 
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FIG. 6. (Color Online). We plot the final system size noo(p) a,s 
a function of the attack size p, when L\ ~ ~ L m = E[L]. 

For the numerical results, we take N = 100, 000, E [L] = 30, 
and consider four tolerance parameters a = 0.2, 0.7, 1.0, 2.0. 
Each data point (represented by a symbol) is the result of 
averaging over 500 independent runs. The lines correspond 
to analytical results computed directly from m- We see a 
perfect agreement between analysis and experiments. In all 
cases, the system breakdowns abruptly through a first order 
transition at p* = — 

^ Ql+l 


is reached, followed by a total breakdown at the critical 


attack size. Numerical results are presented and confirm 
the analytical findings. Last but not least, we prove that 
with mean load fixed, robustness of the power system is 
maximized when the variation among the line loads is 
minimized. In other words, a Dirac delta load distribu¬ 
tion leads to the optimum robustness. 

Our results highlight how different parameters of the 
load distribution and the power line capacity affect the 
robustness of the power grid against failures and attacks. 
To that end, our results can help derive guidelines for 
the robust design of the power grid. We believe that the 
results presented here give very interesting insights into 
the cascade processes in power grids, although through 
a very simplified model of the grid. The obtained results 
can be useful in other fields as well, where equal redis¬ 
tribution of flows is a reasonable assumption. Examples 
include traffic jams, landslides, etc. 

There are many open problems one can consider for 
future work. For instance, the analysis can be extended 
to the case where the tolerance parameter a is not the 
same for all lines, but follows a given probability distri¬ 
bution. It would be interesting to see if the robustness is 
still maximized with a narrow distribution of a. It may 
also be of interest to study robustness against targeted 
attacks rather than random failures. 

ACKNOWLEDGMENTS 

This research was supported in part by National Sci¬ 
ence Foundation through grant CCF ^1422165, and by 
the Department of Electrical and Computer Engineering 
at Carnegie Mellon University. 


[1] M. Rosas-Casals and R. Sole, International Journal of 
Electrical Power & Energy Systems 33, 805 (2011). 

[2] G. Andersson, P. Donalek, R. Farmer, N. Hatziar- 
gyriou, I. Kamwa, P. Kundur, N. Martins, J. Paserba, 
P. Pourbeik, J. Sanchez-Gasca, et al, Power Systems, 
IEEE Transactions on 20, 1922 (2005). 

[3] G. Zhang, Z. Li, B. Zhang, and W. A. Halang, Physica 
A: Statistical Mechanics and its Applications 392, 3273 
(2013). 

[4] Y. Tang, G. Bn, and J. Yi, in Zhongguo Dianji 
Gongcheng Xuebao(Proceedings of the Chinese Society of 
Electrical Engineering), Vol. 32 (Chinese Society for Elec¬ 
trical Engineering, 2012) pp. 167-174. 

[5] S. V. Buldyrev, R. Parshani, G. Paul, H. E. Stanley, and 
S. Havlin, Nature 464, 1025 (2010). 

[6] D. J. Watts, Proceedings of the National Academy of 
Sciences 99, 5766 (2002). 

[7] R. Kinney, P. Crucitti, R. Albert, and V. Latora, The 
European Physical Journal B-Condensed Matter and 
Complex Systems 46, 101 (2005). 

[8] T. D. O’Rourke, BRIDGE-WASHINGTON-NATIONAL 
ACADEMY OF ENGINEERING- 37, 22 (2007). 

[9] I. Dobson, B. A. Carreras, V. E. Lynch, and D. E. New¬ 
man, Chaos: An Interdisciplinary Journal of Nonlinear 
Science 17, 026103 (2007). 


[10] O. Yagan, D. Qian, J. Zhang, and D. Cochran, IEEE 
Transactions on Parallel and Distributed Systems 23, 
1708 (2012), arXiv:1201.2698v2. 

[11] S. Pahwa, C. Scoglio, and A. Scala, Scientific reports 4 
(2014). 

[12] H. Daniels, Proceedings of the Royal Society of London. 
Series A. Mathematical and Physical Sciences 183, 405 
(1945). 

[13] J. V. Andersen, D. Sornette, and K.-t. Leung, 

Phys. Rev. Lett. 78, 2140 (1997) 

[14] T. J. Overbye, X. Cheng, and Y. Sun, in Proceedings, 
37th Hawaii International Conference on System Sci¬ 
ences (2004). 

[15] B. Stott, J. Jardim, and O. Alsac, 

Power Systems, IEEE Transactions on 24, 1290 (2009) 

[16] A. L. Barabasi and L. Albert, Science 286, 509 (1999). 

[17] W. Curtin, Journal of the Mechanics and Physics of 
Solids 41, 217 (1993). 

[18] F. Kun, M. Costa, R. Costa Filho, J. Andrade Jr, 
J. Soares, S. Zapperi, and H. Herrmann, Journal of 
Statistical Mechanics: Theory and Experiment 2007, 
P02003 (2007). 

[19] D. Cohen, P. Lehmann, and D. Or, Water resources 
research 45 (2009). 















12 


[20] S. Pradlian and B. K. Chakrabarti, International Journal [29] 
of Modern Physics B 17, 5565 (2003). 

[21] A. E. Motter and Y.-C. Lai, 

Phys. Rev. E 66, 065102 (2002), [30] 

[22] W.-X. Wang and G. Chen, 

Phys. Rev. E 77, 026101 (2008) 

[23] B. Mirzasoleiman, M. Babaei, M. Jalili, and M. Safari, 

Physical Review E 84, 046114 (2011). [31] 

[24] P. Crucitti, V. Latora, and M. Marchiori, 

Phys. Rev. E 69, 045104 (2004), [32] 

[25] R. da Silveira, Phys. Rev. Lett. 80, 3157 (1998) 

[26] D. Sornette, K.-T. Leung, and J. Andersen, arXiv 

preprint cond-mat/9712313 (1997). [33] 

[27] C. Roy, S. Kundu, and S. Manna, Physical Review E 

91, 032103 (2015). [34] 

[28] V. Rosato, L. Issacharoff, F. Tiriticco, [35] 
S. Meloni, S. Porcellinis, and R. Setola, 


International Journal of Critical Infrastructures 4, 63 (2008) [36] 


The case of non-uniform initial loads in democratic fiber 
bundle model is briefly discussed in m in the context of 
power systems, and some numerical results are provided. 
1. Dobson, J. Chen, J. Thorp, B. A. Carreras, and D. E. 
Newman, in System Sciences, 2002. HICSS. Proceedings 
of the 35th Annual Hawaii International Conference on 
(IEEE, 2002) pp. 10-pp. 

For convenience, we assume that a line also fails when its 
load equals its capacity. 

M. Faloutsos, P. Faloutsos, and C. Faloutsos, in ACM 
SICCOMM Computer Communication Review, Vol. 29 
(ACM, 1999) pp. 251-262. 

A. Clause!, C. R. Shalizi, and M. E. J. Newman, 
SIAM Rev. 51, 661 (2009) 

M. E. J. Newman, Phys. Rev. E 66 (2002). 

E. A. Leicht and R. M. D’Souza, 
arXiv:0907.0894vl [cond-mat.dis-nn] (2009), 0907.0894 
A. Papoulis and S. U. Pillai, Probability, random vari¬ 
ables, and stochastic processes, 4th ed. (Tata McGraw- 
Hill Education, 2002). 


